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ABSTRACT 

Spectrum sensing is a fundamental component in cognitive 
radio. A major challenge in this area is the requirement of a 
high sampling rate in the sensing of a wideband signal. In this 
paper a wideband spectrum sensing model is presented that 
utilizes a sub-Nyquist sampling scheme to bring substantial 
savings in terms of the sampling rate. The correlation matrix 
of a finite number of noisy samples is computed and used by 
a subspace estimator to detect the occupied and vacant chan- 
nels of the spectrum. In contrast with common methods, the 
proposed method does not need the knowledge of signal prop- 
erties that mitigates the uncertainty problem. We evaluate the 
performance of this method by computing the probability of 
detecting signal occupancy in terms of the number of sam- 
ples and the SNR of randomly generated signals. The results 
show a reliable detection even in low SNR and small number 
of samples. 

Index Terms — Wideband spectrum sensing, Sub-Nyquist 
sampling, Cognitive radio, Correlation matrix, Subspace 
methods 

1. INTRODUCTION 

Spectrum sensing is an important function to enable cogni- 
tive radios to detect the underutilized spectrum licensed to the 
primary systems and improve the overall spectrum efficiency 
CI. Some well-known spectrum sensing techniques are en- 
ergy detection, matched filter and cyclostationary feature de- 
tection that have been proposed for narrowband sensing. In 
these methods, based on the signal properties, a decision is 
made to detect presence or absence of a primary user in the 
considered band [ 1 J. 

Future cognitive radios should be capable of scanning a 
wideband of frequencies, in the order of few GHz |2J. In the 
wideband regime, the radio front-end can employ a bank of 
band-pass filters to select a frequency band and then exploit 
the existing techniques for each narrowband, but this method 
requires a large number of RF components . 

Alternatively, in order to identify the locations of vacant 
frequency bands, the entire wideband is modelled as a train of 
consecutive frequency sub-bands HI and the total wideband is 



sampled classically or in a compressed way. After estimation 
of the spectrum from the obtained samples, the conventional 
spectrum sensing methods such as energy detection would 
be applied to detect the signal in each band. Classical sam- 
pling of a wideband signal needs high sampling rate ADCs, 
which have to operate at or above the Nyquist rate. Clearly, 
this is a major implementation challenge. Recent work based 
on compressive sampling has been proposed to overcome the 
problem of high sampling rates in [2], [3]. However, estimat- 
ing the spectrum of a signal from its compressed samples is 
achieved by solving an optimization problem [2 1, which is not 
an easy task. In addition, detection of the signal in each chan- 
nel from its estimated spectrum needs knowledge of the signal 
and noise power that is not reliable because of uncertainty. 

By using the fact that the wireless signals in open- 
spectrum networks are typically sparse in the frequency 
domain, we propose a wideband spectrum sensing method 
that would bring substantial saving in terms of the sampling 
rate. The frequency band of interest is divided into a finite 
number of spectral bands, and the presence and absence of 
the signal in each spectral band is examined by considering 
the correlation matrix of the sampled data. In our method, 
estimation of the signal spectrum is skipped, and we directly 
detect the occupied channels from the sampled data in the 
time domain. Additionally, the problem of noise uncertainty 
is not of concern in this method, since the knowledge of sig- 
nal and noise power is not needed. We evaluate this method 
by computing the probability of detecting signal occupancy 
in terms of the number of samples and signal to noise ratio 
(SNR). The outline of the paper is as follows: The next sec- 
tion states the signal model and problem formulation. Section 
[3] introduces the spectrum sensing method and explains the 
functionality of each block in the model. In Section HI the 
simulation results and performance evaluation are presented 
and finally a conclusion is given in Section|5] 

2. PROBLEM STATEMENT 

The received signal x(t) is assumed to be an analog wideband 
sparse spectrum signal, bandlimited to [0, B max ]. Denote the 
Fourier transform of x(t) by X(f). Depending on the ap- 
plication, the entire frequency band is segmented into L nar- 
rowband channels, each of them with bandwidth B, such that 



B m ax = L x B. It is assumed that the signal bands are un- 
correlated with each other. The channels are indexed from 
to L — 1. Those spectral bands which contain part of the sig- 
nal spectrum are termed active channels, and the remaining 
bands are called vacant channels. Denote the number of such 
active channels by N. The indices of the N active channels 
are collected into a vector 



b = [h, b 2 , ■ ■ . ,6jv] 



(1) 



which is referred to as the active channel set. 

In the considered system, N and b are unknown. How- 
ever, we know the maximum channel occupancy which is de- 
fined as 

N> max 



L 



(2) 



where N max > N is the maximum possible number of occu- 
pied channels. Figure OJdepicts the spectrum of a multiband 
signal at the sensing radio, which contains L = 32 chan- 
nels, each with a bandwidth of B = 10 MHz. The signal 
is present in N = 6 channels, and the active channel set is 
b= [8,16,17,18,29,30]. 

The problem is, given B max , B and !!„„, to find the 
presence or absence of the signal in each spectral band or 
equivalently find the active channel set, b, at a sub-Nyquist 
sample rate. 

3. WIDEBAND SPECTRUM SENSING MODEL 

The proposed model for wideband spectrum sensing is illus- 
trated in Figure |2] The analog received signal at the sens- 
ing cognitive radio is sampled by the multicoset sampler at a 
sample rate lower than the Nyquist rate. The sampling reduc- 
tion ratio is affected by the channel occupancy and multicoset 
sampling parameters. The outputs of the multicoset sam- 
pler are partially shifted using a multirate system, which con- 
tains the interpolation, delaying and downsampling stages. 
Next, the sample correlation matrix is computed from the fi- 
nite number of obtained data. Finally, the correlation matrix 
is investigated to discover the position of the active channels 
by subspace methods. In this section each block of the model 
is described in detail. 

3.1. Multicoset sampler 

The analog wideband signal x(t) is sampled using a mul- 
ticoset sampling scheme introduced in [4|. The multicoset 
sampler provides p data sequences for i — 1, p, given by 



Xi(m) = x[{mL + Ci)/B max ],m e Z, 



(3) 



where {q}, is a randomly selecting p numbers out of the set 

L={o,i,...,£-i}ii,iia. 

The average sample rate of this scheme is f avg — aB max , 
[4 1 where a = (j-) is termed the sub-Nyquist factor. Accord- 
ing to Landau's lower bound J6J, a is lower bounded to the 
maximum channel occupancy, a > Q max . 
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Fig. 1. Spectrum of a wideband signal received at the sensing 
radio with L = 32 total bands and N = 6 active channels. 
The active channel set is b = [8, 16, 17, 18, 29, 30]. 

The number of data sequences, p, should be chosen 
greater than the maximum number of active channels N max , 
to satisfy the Landau's lower bound and provide enough 
equations to find the unknown parameters. 

3.2. Sample correlation matrix 

The main purpose of this section is to relate the problem of 
spectrum sensing with the problem of parameter estimation. 
Towards this goal, the correlation matrix of a special config- 
uration of sampled data is computed. In order to achieve this, 
the following configurations are applied on the sampled data. 

First, each Xi (m) sequence is over-sampled by a factor L, 
such that 

x,i(j), n = mL, m € Z 



0, 



otherwise 



and then it is filtered to obtain x\ li [n] = x Ui [n] * h[n], where 
h[n] is the interpolation filter with the frequency response of 



H(f) 



1, f€[0,B] 
0, otherwise. 



(4) 



Next, the output filtered sequence is delayed with q samples 
such that 

x Ci [n] = Xhi [n - Ci] . (5) 
Let us define y(/) as the known vector of observations 



y(f)=[X 1 (f),X 2 (f),...,X p (f)\ 



(6) 



where the superscript T denotes the transpose, and Xi (/) is 
the DFT of the sequence x Ci [n]. Also, x(/), the unknown 
vector of the signal spectrum parameters is defined as 



x(/) = 



X{f + b x By 
X(f + b 2 B) 

X(f + b N B) 



fe[o,B] 



(7) 



where X(f + biB), f S [0, B], are the frequency elements of 
the signal in the active band indexed by 6j. 
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Fig. 2. Proposed wideband spectrum sensing model. 



After applying Fourier transform on both sides of (0 and 
expressing the result in matrix form, the data model in the 
frequency domain is given by [5 ] 

y(/)=A(b)x(/)+n(/), fe[0,B] (8) 

where A(b) e cp xW j s the modulation matrix given by [5] 

'j2ircib k 



A(b)(i,fc) = Bexp 



L 



(9) 



and n(/) is the frequency representation of the noise. For 
simplicity, we assume that n(/) is a Gaussian complex noise 
with distribution of A/"(0,<t 2 I), which is also uncorrected 
with the signal. 

The model in ^ is a classical signal model that relates 
the observation vector, y(/), with the unknown signal spec- 
trum vector, x(/), via the modulation matrix A(b). Note that 
the unknown signal parameter b is also the active channel set 
that is desired in the problem of spectrum sensing. Therefore, 
with this configuration, the problem of wideband spectrum 
sensing is turned into the problem of finding the model pa- 
rameter b with minimum length N, subject to the data model 
(H}. This is a combined detection-estimation problem, where 
we want to estimate the number as well as the parameters of 
the signals. Due to the computational complexity involved, 
the problem is solved in two steps: first the number of param- 
eters is detected, and then, with an estimate of the number of 
parameters N at hand, the parameters of the signal are esti- 
mated. An approach based on the correlation matrix of the 
observations is employed here for solving this problem. 

The correlation matrix of observation vector is defined as 



R = £[y(/)y*(/)] = A(b)PA*(b) 



Vl 



where ()* denotes the Hermitian transpose, and 

P = £[x(/)x*(/)] 



(10) 



(ID 



is the correlation matrix of the signal vector [7 1. 

Since the distribution of the signal is unknown the real 
correlation matrix R cannot be achieved. Hence, we estimate 
R from the integration of [y (/)y* (/)] over the interval [0, B\. 
However, from Parseval's identity it can be computed directly 
in the time domain from the sequences x Ci [n] at the sample 
rate of B max . Since each x Ci [n] sequence is the output of a 



narrowband filter, the reduced bandwidth output signal can be 
easily accommodated within a lower output sample rate. This 
means that the computations do not need to be performed at 
the high sample rate, B max . Thus, the sequences are down- 
sampled by L, the reduced bandwidth factor, such that 



Xdi (jn) — x Ci [mL] 



(12) 



The total process of oversampling, filtering, delaying and 
downsampling from Xi(m) to x^m) is viewed as a frac- 
tional shifting of the sequence Xi(n). Defining the snapshot 
vector Xd(m) as 



x d (m) 



x dl (m) 
x d2 (m) 

Xd p (m) 



the p x p sample correlation matrix from M samples of the 
partially shifted sequence is computed from the formula Q 



1 M 

R= *d( m )x d (m). 



Under suitable assumptions R > R when M — > oo. 



(13) 



3.3. Subspace Analysis 

The signal in each spectral band is assumed to be uncorrected 
with the other bands. Hence, the correlation matrix of the 
signal vector, P, is always full-rank and the class of subspace 
methods is applicable which is described in this section. 

3.3.1. Estimating the Number of Active channels 

In a geometrical view of ( [Tol l, any orthogonal vector to A(b) 
is an eigenvector of R with corresponding eigenvalue a 2 . The 
remaining eigenvectors are all in the range space of A(b), 
and therefore are termed signal eigenvectors. The eigen- 
decomposition of R is partitioned into a signal and a noise 
subspace as Q 



R — E.,A.,E* 



E n A„E* 



(14) 



where A s and A„ are diagonal matrices of signal and noise 
eigenvalues respectively, and E s and E„ are the matrices of 
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Fig. 3. Ordered eigenvalues of a sample correlation matrix 
with M = 61. There are six eigenvalues greater than the 
threshold level. 



the corresponding eigenvectors. The signal eigenvectors in 
E s span the range space of A(b) , which is termed the signal 
subspace. For the noise eigenvector we have instead, E„ _L 
A(b) [7 |. So, if the dimension of noise subspace is obtained, 
the orthogonality property can be exploited to find the signal 
parameters. 

If we denote the ordered eigenvalues of R by 



A, > A 2 > 



then it follows that in case of large enough number of samples, 
R > R and the smallest (p — N) eigenvalues of R are all 
equal to a 2 0. Thus, for large M, the dimension of the signal 
vector can be determined from the multiplicity of the smallest 
eigenvalues of R. 

In realistic scenarios, the number of samples depends on 
the duration of the sensing period, which should be as short as 
possible, especially in case of time-varying channels. Hence 
the sample matrix R ^ R. In the spectrum of ordered eigen- 
values, the signal eigenvalues are still identified as the N 
largest ones. But, the noise eigenvalues are no longer equal to 
each other (9) and the separation between the signal and noise 
eigenvalues needs a threshold level. Depending on the noise 
power and the number of samples, the threshold level should 
be changed. Figure [3] depicts a typical set of ordered eigen- 
values that are separated by a horizontal line as the threshold 
level. As one can see there are six eigenvalues greater than the 
threshold level. That specifies the number of active channels. 

To avoid the threshold setting for different tests, we ap- 
plied information theoretic criteria for model order selection 
such as minimum description length (MDL). The number of 
active channels using the MDL criterion for < r < N max 



Fig. 4. Frequency representation, X(f), and the correspond- 
ing Pmu values of a typical wideband signal with L = 32 
channels. The position of six significant values specify the 
occupied channels b = [8, 16, 17, 18, 29, 30]. 



is given by MM 

N = argmin — M (p — r) log 



g( r ) 

o(r) 



-\r{2p- 



ilogM (15) 



where M is the number of samples in each sequence, g{r) 
and a(r) are the geometric and arithmetic mean of the (p — r) 
eigenvalues of the correlation matrix respectively. 

The probability of correct detection of N, namely, Pr(N = 
N), depends on the number of samples (M), SNR and sub- 
Nyquist factor (a). Since evaluation of this detection al- 
gorithm in the Neyman-Pearson sense of finding the most 
powerful test that does not exceed a threshold probability 
of false alarm is computationally very expensive J8), we 
will evaluate the performance of this method numerically in 
SectionS] 

3.3.2. Active channel set recovery 

After estimating the number of active channels, the (p — N) 
smallest eigenvalues are specified as the noise eigenvalues. 
Denote the corresponding noise eigenvectors with E„ as a 
p x (p ~ N) matrix. The location of the active channels can 
be recovered according to a MUSIC-Like algorithm as 1 10 1 

1 



P M u(k) 



|a fe E r , 



< k < L - 1 



(16) 



where ||.|| denotes the 2-norm, k is the channel index and afe 
is a column of A(b), given by 



afe 





e £ , e £ 



(17) 



The algorithm generates L values corresponding to the L 
channels. If k is the index of an active channel, PMu(k) is 
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Fig. 5. Pr(N = 6), probability of detecting number of ac- 
tive channels versus M and SNR for the simulated wideband 
system. 

significant in that point, otherwise it will be a small value. The 
estimated active channel set, b, is determined by selecting the 
position of these significant values. Figure |4]depicts the fre- 
quency representation of the received signal, X(f), and the 
corresponding Pmu values of a typical wideband system. As 
illustrated in the figure, several significant values correspond 
to active channels are appeared where their locations specify 
the estimated active channel set. The other channels are inter- 
preted as the vacant channels and can be used by the cognitive 
system to transmit. 

4. NUMERICAL RESULTS 

In this section, we illustrate the performance of our method 
using Monte Carlo simulations. 

The received signal at the cognitive radio sensing is gen- 
erated from the model 



JV 

£■ 

i=l 



n[n) * h[n])exp(j2irfin/B max ) + w[n] 



where * shows a convolution between ri[n] ~ A/"(0, of) and 
h[n], the low pass filter defined in ©. The output of the con- 
volution is placed at the carrier frequency and corrupted 
by w[n] ~ Af(0, 1), the additive white Gaussian noise. 

The wideband of interest is in the range of [0, 320] MHz, 
containing 32 channels of equal bandwidth of B = 10MHz. 
The signal variance is chosen such that the received SNR of 
all active channels are equal. Figure [T]depicts the spectrum of 
the signal model with N = 6 active bands located at different 
unknown carriers. Given B max = 320MHz, Vl max = 0.25 
and B = 10MHz, it is desired to find the positions of occu- 
pied and vacant channels at a sub-Nyquist sampling rate. 



Fig. 6. Detection performance of the proposed model versus 
M and SNR for the simulated wideband system.. 



A multicoset sampler with parameters L = 32, p = 10 
is used to sample the signal at the average sample rate of 
favg = 100 MHz, which is a « 0.3 of the Nyquist rate. 
Ten Ci numbers are selected randomly out of the set L. 

We first investigate the performance of estimating N. Us- 
ing 1000 Monte Carlo simulations, and various values of M 
and SNR, we obtain an estimate of the number of active chan- 
nels from the eigenvalues of the sample correlation matrix us- 
ing <fl5V For a choice of SNR and different values of M, 
we compute the empirical probability of detecting six active 
channels. As M and SNR are increased, we expect the esti- 
mator to detect the correct value with high probability. 

Figure [5] depicts the computed Pr(N = 6) for different 
number of samples and SNR. It is seen at SNR=ldB after 
M > 41, the estimator is able to detect the correct value with 
high probability. However, in lower SNR, more samples are 
needed to achieve a high probability. 

Next, we evaluate the detection performance of the pro- 
posed method by computing the probability of detecting the 
signal occupancy as 



1 N 

p d = ^E Pr ^ eS i^ eb ) 



and the false alarm probability as 



(18) 



Pf 



1 



L-N 



L-N 

E 

i=l 



Pr{b\ e \>\b\ e b c ) 



(19) 



where b c = L b is the complement set of b. 

The P^ and Pe are computed from ( fT8l and ( fT9] i for dif- 
ferent values of M and SNR and the results are shown in Fig- 
ure [6] The results exhibit outstanding detection performance 
even in low SNR and small M. It is seen that at SNR=ldB, 



after M > 41, the proposed model detects the occupied chan- 
nels with probability close to one. The value of Pf decreases 
dramatically with increasing M, such that after M > 21 it 
becomes almost zero for all SNR values. Comparing Fig- 
ure|5]and [6]reveals that with correct estimation of N, the per- 
fect detection of channel occupancy is possible. In the other 
words, the accuracy of the proposed model only depends on 
the MDL estimator and would be improved using other tech- 
niques such as the method in (9) . 

5. CONCLUSION 

A method of wideband spectrum sensing for cognitive radio 
is proposed to mitigate the limitations of high sampling rate, 
high complexity and noise uncertainty. The proposed tech- 
nique utilizes a multicoset sampling scheme that can use ar- 
bitrarily low sampling rates close to the channel occupancy. 
With low spectrum utilization assumption, this would bring 
substantial savings in terms of the sampling rate. The coset 
samples are fractional shifted and used to compute the corre- 
lation matrix of the signal. The computation cost of this step 
is linear in the amount of data. The problem of spectrum sens- 
ing is turned into the problem of parameter estimation and 
then is solved by subspace methods that today are standard 
tools in signal processing. We evaluate the detection perfor- 
mance of this method for a typical case. The results show 
that even in low SNR with taking enough number of samples 
a perfect detection is possible. For a typical wideband sys- 
tem with il max = 0.25 and SNR=ldB by taking M = 31 
samples, at a ~ 0.3 of the Nyquist rate, P^ = 0.99 and 
Pf = 10 -3 are achieved. 
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